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Abstract — We study a problem where wireless service providers 
compete for heterogenous wireless users. The users differ in 
their utility functions as well as in the perceived quality of 
service of individual providers. We model the interaction of an 
arbitrary number of providers and users as a two-stage multi- 
leader-follower game. We prove existence and uniqueness of 
the subgame perfect Nash equilibrium for a generic channel 
model and a wide class of users' utility functions. We show 
that the competition of resource providers leads to a globally 
optimal outcome under mild technical conditions. Most users 
will purchase the resource from only one provider at the unique 
subgame perfect equilibrium. The number of users who connect 
to multiple providers at the equilibrium is always smaller than the 
number of providers. We also present a decentralized algorithm 
that globally converges to the unique system equilibrium with 
only local information under mild conditions on the update rates. 

Index Terms — Game theory, provider competition, wireless 
network, pricing 



I. Introduction 

Due to the deregulation of the telecommunication industry, 
future wireless users are likely to freely choose a provider 
(or providers) offering the best tradeoff of parameters. This 
is already happening with some public Wi-Fi connections, 
where users can connect to wireless access points of their 
choice, with usage-based payments and no contracts. Despite 
the common presence of a free public wi-fi network, some 
users may still choose more expensive providers who offer 
better quality of service. Another example of the deregulation 
trend is the analog television (UHF) spectrum which was 
recently open for unregulated use in the US [ 1 1. 

In this work, we consider a situation where wireless service 
providers compete to sell a limited amount of wireless re- 
sources (e.g., frequency bands, time slots, transmission power) 
to users who are free to choose their provider(s). We in- 
vestigate how providers set prices for the resource, and how 
users choose the amount of resource they purchase and from 
which providers. The focus of our study is to characterize 
the outcome of this interaction. We consider the general 
case where different users have different utility functions and 
experience different channel conditions to different service 
providers. 

We model the user-provider interaction as a multi-leader- 
follower game (see p}, pi). The providers announce the 
wireless resource prices in the first stage, and the users 
announce their demand for the resource in the second stage. 



A user may purchase resource from a provider with a poor 
channel if the price of the resource is low, or conversely from 
an expensive one if the channel is strong. The providers select 
their prices to maximize their revenues, keeping in mind the 
impact of their prices on the demand of the users. 
The contributions of our work are as follows: 

• General Heterogeneous Wireless Network Model: We 
study a general network model that captures the hetero- 
geneity of wireless users and service providers. The users 
have different utility functions, the providers have differ- 
ent resource constraints, the channel gains between users 
and providers are independent and arbitrarily disttibuted, 
and the numbers of users and providers can be arbitrary. 

• Unique Socially Optimal Allocation: We first study the 
problem of maximizing social welfare under a fairly 
general utility function model. We show that when chan- 
nel parameters are randomly drawn from continuous 
distributions, there exists a unique optimal solution to 
the problem with probability 1, despite the non-strict 
convexity of the optimization problem. 

• Existence, Uniqueness, and Zero Efficiency Loss of Equi- 
librium: We further prove existence and uniqueness of 
the subgame perfect Nash equilibrium in the two-stage 
game, under an easily verifiable sufficient condition on 
the users' utility functions. Moreover, we show that the 
unique equilibrium maximizes the social welfare, despite 
the selfish nature of the providers and users. 

• Primal-Dual Algorithm converging to Equilibrium: We 
provide a decentralized algorithm that results in an equi- 
librium of the provider competition game. The partici- 
pants only need local information during the execution of 
this algorithm. Providers only need to know the demand 
of the users, while users only need to consider the prices 
given by the providers. 

We begin by describing the provider competition model 
in Section III] In Section III we discuss the socially optimal 
resource allocation and in Section IV we analyze the provider 
competition game. We define the primal-dual update algorithm 
and prove its convergence in Section W\ We provide numerical 
results and dicussion in Section IVll We discuss the related 
work in Section IVIII and conclude in Section IVIIII 

II. Problem Formulation 

We consider a set J = {1, . . . , J} of service providers and 
a set I = {!,...,/} of users. Provider j 6 J maximizes its 



revenue by selling up to Qj amount of resource to the users. 
A user i el maximizes its payoff by purchasing resources 
from one or more providers. The communication can be both 
downlink or uplink, as long as users do not interfere with each 
other by using orthogonal resources. We model the interaction 
as a multi-leader-follower game (see pj, (3)), where providers 
are the leaders and users are the followers. We called this 
game the provider competition game. We assume that the users 
play the game during the coherence time of their channels. 
This can be reasonable for a quasi-static network environment 
(e.g., users are laptops or smart phones in offices or airports). 
This means that the channel gains remain roughly constant 
and can be made known to all parties. For example, each 
provider collects its channel condition information to each 
user, and then broadcasts this information to all users. This 
assumption will be relaxed in Section [Vl where we consider 
a decentralized algorithm that results in the same outcome as 
the multi-leader-follower game. 

A. Provider competition game 

The provider competition game consists of two stages. In 
the first stage, providers announce prices p — [pi ■ ■ ■ p.j], 
where pj is the unit resource price charged by provider j. In 
the second stage, each user i £ I chooses a demand vector 
Qi = [qh ■ ■ ■ Qij], where qij is the demand to provider j. We 
denote by q = [q-i • ■ ■ qj\ the demand vector of all users. 

In the second stage where prices p are given, the goal of 
user i is to choose q i to maximize its payoff, which is utility 
minus payment: 



Vi(qi,p) 




(1) 



where c,j is the channel quality offset for the channel between 
user i and the base station of provider j (see Example [T] and 
Assumption E), and m is an increasing and concave utility 
function. In the first stage, a provider j chooses price pj 
to maximize its revenue pj Yli=i Qij subject to the resource 
constraint 2~2i=i Qij — Qj> by taking into account the demand 
of the users in the second stage. We consider linear pricing 
with no price discrimination across the users. 

Under this model, a user is allowed to purchase from 
several providers at the same time. For this to be feasible, 
a user's device might need to have several wireless interfaces. 
Mathematically, the solution of this model gives an upper 
bound on best performance of any situation where users are 
constrained to purchase from one provider alone. Interestingly, 
our results show that for most users, i.e. no less than I — J, 
the optimal strategy is to choose exactly one provider. 

Next we give a concrete example of how our model is 
mapped into a physical wireless system. 

Example 1: Consider wireless providers operating on or- 
thogonal frequency bands Wj, j G J. Let q^ be be the 
fraction of time that user i is allowed to transmit exclusively 
on the frequency band of provider j, with the constraint 
Tliei- Q^ = 1' i ^ 3 ' Furthermore, assume that each 
user has a peak power constraint Pj. We can then define 



cij — Wj log(l H — i U) ), where hij is the channel gain and 

of, is the Gaussian noise variance for the channel between 
user i and network j. In this case, a user's payoff is the 
difference between its utility function (in terms of total rate) 
and payments, v t = u»(£)j=i q%jCij) - J2j=iPjQtj- 

Although the c^ channel quality offset factor represents 
channel capacity in Example [T] it can be any increasing 
function of the channel strength depending on the specific 
application scenario. 

Finally, we remark that the problem shares certain similarity 
with the multi-path routing problem in a generalized network 
flow setting, where each source corresponds to a user and each 
link corresponds to a provider. The key difference is that in our 
model the user-provider connections have different weights, 
which is not the case for the multipath routing problem. 

B. Model assumptions 

We make the following assumptions throughout this paper: 

Assumption 1: (Utility functions): For every user i £ X, 
Ui(x) is differentiable, increasing, and strictly concave in x. 
This is a standard way to model elastic data applications in 
network literature (see, e.g., [4|). 

Assumption 2: (Channel quality offsets and channel gains): 
Channel quality offsets c^ are drawn independently from 
continuous, possibly different utility distributions. In particular 
Pr(cij = Cki) = for any i, k € X, j,l e J. The channel 
quality offset accounts for the effect that buying the same 
amount of resource from different providers typically has 
different effects on a user's quality of service. As Example 
fl] shows, channel quality offset Cy may be a function of 
the channel gain hij between user i and provider j. In this 
case the assumption is fulfilled if channel gains are drawn 
from independent continuous probability distribution (e.g., 
Ralyleigh, Rician, distance-based path-loss model). 

Assumption 3: (Atomic and price-taking users): The de- 
mand for an atomic user is not infinitely small and can have 
an impact on providers' prices. Precise characterization of this 
impact is one of the focuses of this paper. On the other hand, 
users are price-takers by the assumption of the two-stage game, 
and do not strategically influence prices. 

To analyze the properties of the provider competition game, 



in Section III we study a related socially optimal resource 
allocation problem and show the uniqueness of its solution in 
terms of users' demands. Then, in Section \W\ we come back 
to the provider competition game. We show that the unique 
socially optimal solution solution corresponds to the unique 
equilibrium of the provider competition game, in which case 
the selfish and strategic behavior of providers and users leads 
to zero efficiency loss. 

III. Social Optimum 

A. Social welfare maximization 

In this section we consider a social welfare problem, which 
aims at maximizing the sum of payoffs of all participants, 
(users and providers). The social welfare problem is equivalent 



to maximizing the sum of users' utility functions since the 
payments between users and providers cancel out. We show 
the uniqueness of its solution in terms of users' demands. For 
clarity of exposition, we define the following notation. 

Definition 1: (Effective resource) Let x = [x\ • ■ ■ Xj] be 
the vector of effective resources, where Xi(q i ) = X}i=i 1ij c ij 
is a function of user i's demand q i = [qn . . . qtj]. 

The social welfare optimization problem (SWO) is: 



SWO : max u(x) = 2_] "i (a 



1=1 



subject to y, Qij c ij = Xi i (zl 
i 

Z^Qij = Qj, 3 e J 

8=1 

over qij , Xi > Vi G I, j e J . 



(2) 

(3) 

(4) 
(5) 



We expressed the SWO in terms of two different variables: 
effective resource vector x and demand vector q, even though 
the problem can be expressed entirely in terms of q. In 
particular, a vector q uniquely determines a vector x through 
equations d3|, i.e. we can write x as x(q). With some abuse 
of notation we will write u(q) when we mean u(x(q)). 

Lemma 1: The social welfare optimization problem SWO 
has a unique optimal solution x* . 

Proof: Since Ui(xi) is strictly concave in Xi, then u(x) = 
2~2i=i u i( x i) i s strictly concave in x. The feasible region 
defined by constraints (|3]l-((5| is convex. Hence, u(x) has a 
unique optimal solution x* subject to constraints (pl-pl. ■ 



B. Uniqueness of the socially optimal demand vector q* 

Even though Ui()'s are strictly concave in Xi, they are 
not strictly concave in the demand vector qj. Hence, SWO 
is non-strictly concave in q. It is well-known that a non- 
strictly concave maximization problem might have several 
different global optimal optimizers (several different demand 
vectors q in our case) (see e.g. |5J,||6)). In particular, one can 
choose Cm's, Qj's, and Uj(-)'s in such a way that a demand 
maximizing vector q* of SWO is not unique. However, we 
can show that such cases arise with zero probability whenever 
channel offsets factors Cy 's are independent random variables 
drawn from continuous distributions (see Assumption [2]). 

In the remainder of this section, we show that SWO has 
a unique maximizing demand vector with probability 1. We 
begin by proving Lemma [2] which is an intermediate result 
stating that any two maximizing demand vectors of SWO must 
have different non-zero components. We then observe that any 
convex combination of two maximizing demand vectors is also 
a maximizing demand vector. Finally, we show that all convex 
combinations of maximizing demand vectors have the same 
non-zero components, which is a contradiction with Lemma 
E] This proves the main result of this section (Theorem [3]). 



To make our argument precise, we first define the support 
set of a demand vector q t as follows. 

Definition 2: (Support set): The support set Ji{qi) of a 
demand vector q { contains the indices of its non-zero entries: 

Ji(Qi) = {J£J-- Qij > 0}. 

Given a demand vector q, the ordered collection of support 
sets Ji,J2,---,Ji is denoted by {Ji}{ =l . 
The support set contains providers that user i has strictly 
positive demand from. 

Lemma 2: Let q* be an optimal solution of SWO (a max- 
imizing demand vector) and {J7i}f =1 be the corresponding 
collection of support sets. Then, q* is almost surel>[J the 
unique maximizing demand vector corresponding to { Jj}f =1 . 
Proof: For a maximizing demand vector q* , equations 
(pl-pl hold, with J2j=i itj = x *i- To P rove the lemma, we 
will uniquely construct q* from x* and {Ji\\ =1 . 

We can divide the users into two categories. The decided 
users purchase from only one provider (\Ji\ — 1), and the 
undecided users from several (\Ji\ > 1). It is also possible 
that some users have zero demand to all providers, but without 
loss of generality, we treat such users as decided. Recall that 



for all users we have x* 



E i= i * 



ij^j 



For a decided user 



i who purchases only from provider j, this reduces to x* = 
q*jCij, and the corresponding unique demand vector is q* = 
[0- • -0^0- ■ -01. 

For undecided users, finding the unique q* is less straight- 
forward as there is more than one q { such that J2 6 f. Qij c ij — 
x*. To show that the demand of undecided users is unique, 
we construct the bipartite graph representation (BGR) Q of 
the undecided users' support sets as follows. We represent 
undecided users by circles, and providers of undecided users 
as squares. We place an edge (i,j) between a provider node 
j and a user node i if j e Ji . 

We give an example of a BGR in Fig. Ill where J\ — 
{a, b, c}, J 2 = {b, d}, J 3 = {d, e, /} and J 4 = {6, g}. 

Pi <llb C ±b + qtgCig 
4~" 




Fig. 1: Bipartite graph representation 

The BGR has the following propertied] (see Fig. Ill: 

1) The sum of effective resource on all the edges con- 
nected to user i is the optimal effective resource x* = 
2~2j e j. Qij c ij — Pi- Borrowing from coding theory and 
with some abuse of terminology, we call Pi the check- 
sum of user node i, 

2) The sum of all edges connected to provider node j 
equals to the difference between the supply Qj and the 

'This holds on the probability space defined by the distributions of Cij's. 
2 Fig. 1 shows a connected graph, but this need not be the case. 



demand from decided users who connect to provider j: 



L^ivAi. 



Sj. We call Sj 



3) 



:(i,j)eg Qij ~ "] 2si:(i,j)<£g lij 

the check-sum of provider node j. 

With probability 1, the BGR does not contain any loops. 



This is proved in Appendix A-l 



As it is the case in Fig. [Tl the number of undecided users 
is smaller than the number of providers. This is a direct 
consequence of Property BJ, and will be proved later. 

We can use the BGR to uniquely determine the demands of 
undecided users. Here we use Fig. [T] as an illustrative example, 
the formal description is given in Appendix A-2 We call this 



procedure the BGR decoding algorithm. Consider leaf node (a 
node with only one edge) g and edge q\ . The BGR implies 
that user 4 is the only undecided customer of provider g. Since 
the demands of all decided users have been determined, then 
we know that q\ g — S g — Q g — J2i=£4 itg- We can then remove 
edge q\ and node g from the BGR, and update the check-sum 
value of node 4 to P4 = x\ — q\ g C4 g . Now consider node 4 and 
edge q\ b . Since edge g| fe is now the only edge connecting with 
user node 4, we have q^c^ — P4 and hence q\ b = P^jc^. 
Next we can consider node a, e, or /, and so on. Property 
[3T> is crucial in this procedure since it guarantees that we can 
always find a leaf node in the reduced graph. 

In each step of the algorithm, we determine the unique value 
of q*A associated with the edge of one leaf. This value is 
independent of the order in which we pick the leaf nodes, 



as seen from Appendix A-2 So, we can construct unique 



demand vector q* for each undecided user i. Together with the 
unique demand vectors of the decided users, we have found the 
unique maximizing demand vector x* of SWO with support 

sets {Ji)l=i- ■ 

Theorem 3: The social welfare optimization problem SWO 
has a unique maximizing solution q* with probability 1. 

here 



Proof: The detailed proof is in Appendix A-3 



we provide an outline. Assume there exist two maximizing 
demand vectors of SWO which, by Lemma |2j have different 
supports sets. The support set of a non-trivial convex combi- 
nation of any two non-negative vectors is the union of support 
sets of these two vectors. Hence, all convex combinations 
of two maximizing demand vectors of SWO, which are also 
maximizing demand vectors, have the same support. This is a 
contradiction to Lemma [2] ■ 

Given an optimal demand vector q* of the SWO problem, 
there exists a unique corresponding Lagrangue multiplier vec- 
tor p*, associated with the resource constraints of J providers 
|7). Next, we show that (q*,p*) is the unique equilibrium of 
the provider competition game defined in Section [TTJ 

IV. Analysis of the two-stage game 

In this section we show that there exists a unique equi- 
librium (defined more precisely shortly) of the multi-leader- 
follower provider competition game. In particular, this equi- 
librium corresponds to the unique social optimal solution of 
SWO and the associated Lagrange multipliers. The idea is to 
interpret the Lagrange multipliers as the prices announced by 
the providers. Moreover, we show that there are at most J — 1 
undecided users at this equilibrium. 



First, we define the equilibrium concept (SJ: 
Definition 3: (Subgame perfect equilibrium (SPE)) A price 
demand tuple (p* ,q*(p*)) is a subgame perfect equilibrium 
for the provider competition game if no player has an incentive 
to deviate unilaterally at any stage of the game. In particular, 
each user i £ 1 maximizes its payoff given prices p* . 
Each provider j € J maximizes its revenue given other 
providers' prices p*_ J = (pi,--- ,p*_ 1; ;p* +1 , ' ' "P*j) and me 
users' demand q*{p*)- 

We will compute the equilibrium concept using backward 
induction. In Stage II, we will compute the best response of the 
users q* (p) as a function of any given price vector p. Then 
in Stage I, we will compute the equilibrium prices p*. For 
equilibrium prices p* , the best response of the users q*{p*) 
is uniquely determined via BGR decoding. 

A. Equilibrium strategy of the users in Stage II 

Consider users facing prices p in the second stage. Each 
user solves a user payoff maximization (UPM) problem: 

(J \ J 

UPM : rnaxn, = max«, /,9ij c ij — / PjQij (6) 



\.i =1 



,/=i 



Lemma 4: For each user i e I, there exists a unique 
nonnegative value x*, such that X) 7 =i c ijQij — x \ f° r every 
maximizer q i of the UPM problem. Furthermore, for any j 



such that 



>0, 



p, 



rnintej 



Pk 

£ik ' 



Proof is given in Appendix B-l We remark that, together, 
the unique x*'s from Lemma [4] form a vector x* , which is 
equal to x* , the SWO maximizer from Lemma [T] 

Definition 4: (Preference set) For any price vector p, user 
i's preference set Ji(p) includes each provider j E J with 



El 



mm keJ 



In light of Lemma HI J. L is the set of providers from which 
user i might request a strictly positive amount of resource. 
Users can again be partitioned to decided and undecided based 
on the cardinality of their preference sets, analogous to the 



distinction made in Section III The preference set of a decided 
user i contains a singleton, and there is a unique vector q t 
that maximizes his payoff. By contrast, the preference set of 
an undecided user i contains more than one provider, and any 
choice of q i > such that x* — J2jeJ- 9*i c y maximizes his 
payoff. 

There is a close relationship between the support sets from 



Section III and preference sets defined here. Facing prices p, a 
user i may request positive resource only from providers who 
are in his preference set Ji. By definition, he actually requests 
positive resource from providers who are in his support set Ji. 
So the support set of a user is a subset of his preference set: 
Jiiqip)) C Ji{p)- We can construct a BGR based on the 
preference sets and show that this BGR also has no loops 



with probability 1, following a similar proof (Appendix A-l 1. 

Suppose that the optimal Lagrange multipliers p 
Section 



III 



from 
users have 



are announced as prices. Since all 

access to complete network information, each of them can 

calculate all users' preference sets, and can construct the 



corresponding BGR. Undecided users can now uniquely deter- 
mine their demand vector by independently running the same 
BGR decoding algorithm. The demand found through BGR 
decoding is unique as all demand vectors are considered at 
one time and equality of supply and demand is taken into 
account. We note that the demand found in this way is only 
one of an undecided user's infinitely many best responses 
under prices p* . However, only the demands given by the BGR 
decoding algorithm will balance the supply and demand for 
each provider at the optimal price p* . We will later show that 
this is the only subgame perfect equilibrium of the provider 
competition game. 

B. Equilibrium strategy of the providers in Stage I 

The optimal choice of prices for the providers depends on 
how the users' demand changes with respect to the price, 
which further depends on the users' utility functions. The 
quantity that indicates how a user's demand changes with 
respect to the price is the coefficient of relative risk aversion 
|9| of utility function u t , i.e. k RRA = — x ^,l ■ We focus on 
a class of utility functions characterized in Assumption HI 

Assumption 4: For each user i <E I, the coefficient of 
relative risk aversion of its utility function is less than 1 . 

Assumption HI is satisfied by some commonly used utility 
functions, such as log(l + x) and the a— fair utility functions 
\_ a , for a € (0, 1) j 10 . Under Assumption 4 a monopolistic 
provider will sell all of its resource Qj to maximize its 
revenue. Intuitively, when a provider lowers its price, the 
demand of the users increases significantly enough that the 
change in revenue of the provider is positive. This encourages 
the provider to lower the price further such that eventually total 
demand equals total supply. In the case of multiple providers, 
As sumption HJ also ensures that all providers are willing to sell 
all their resources to maximize their revenues. 

We call the prices that achieve equality of demand and 
supply market clearing prices. 

Theorem 5: Under Assumption HI the unique socially opti- 
mal demand vector q* and the associated Lagrangian multi- 
plier vector p* of the SWO problem constitute the unique sub- 
game perfect equilibrium of the provider competition game. 



The proof is given in Appendix B-2 It is interesting to 
see that the competition of providers does not reduce social 
efficiency. This is not a simple consequence of the strict 
concavity of the users' utility functions; it is also related to 
the elasticity of users' demands. Assumption HI ensures that 
the demands are elastic enough such that a small decrease in 
price leads to significant increase in demand and thus a net 
increase in revenue. 

Under the optimal prices p* announced by the providers in 
the first stage, the users in the second stage will determine 
the unique demand vector q* using BGR decoding. On the 
other hand, if the providers charge prices other than p*, no 
best-response from the users will make the demand equals to 
the supply, which is a necessary condition for an equilibrium. 

In light of Theorem BJ we will refer to the unique subgame 
perfect equilibrium (p* , q*) of the provider competition game 
as the equlibrium. 



C. The number of undecided users 

Since the presence of undecided users makes the analysis 
challenging, it is interesting to understand how many unde- 
cided users there can be in a given game. It turns out that 
such number is upperbounded by the number of providers J 
in the network. 

Lemma 6: Under any given price vector p in the first stage, 
the number of undecided users in the second stage is strictly 
less than J. 

The proof is given in Appendix |B-3| The main idea is that 
if the number of undecided user nodes in a BGR is not smaller 
than the number of provider nodes, then there exists a loop 
in the BGR. This, however, occurs with zero probability, as 
shown in Section IHTl 

V. Primal-Dual Algorithm 

The previous analysis of the subgame perfect equilibrium 
has assumed that every player (provider or user) knows the 
complete information of the system. This may not be true 
in practice. In this section we present a distributed primal- 
dual algorithm where providers and users only know local 
information and make local decisions in an iterative fashion. 
We show that such algorithm globally converges to the unique 
equilibrium discussed in Theorem BJ under mild conditions on 
the updating rates. 

The key proof idea is to show that the primal-dual algorithm 
converges to a set containing the optimal solution of SWO. We 
can further show that this set contains only the unique optimal 
solution in most cases, regardless of the values of the updating 
rates. 

We first present the algorithm, and then the proof of its 
convergence. 



A. Primal-dual algorithm 

In this section we will consider a continuous-time algorithm, 
where all the variables are functions of time. For compactness 
of exposition, we will sometimes write qij and pj when we 
mean (&,-(£) and Pj(t), respectively. Their time derivatives -^ 



and 



■'!': 



at 



~Uf will often be denoted by q^ and pj. We denote by q* 
and p* the unique maximizer of SWO and the corresponding 
Lagrange multiplier vector, respectively. As shown in Theorem 
E\ (p* ,q*) is also the unique subgame perfect equilibrium of 
the provider competition game. These values are constant. 

To simplify the notation, we denote by fij(t) or simply 
fij the marginal utility of user i with respect to q^ when his 
demand vector is qi(t): 



Jij — 



dui(qi) 



dui{x) 



dx 



=^=E/ = i9i. 



(7) 



We will use /*■ to denote the value of fij(t) evaluated at q*, 
the maximizing demand vector of user i. So, /*■ is a constant 
that is equal to a user's equilibrium marginal utility as opposed 
to fij(t) which indicates marginal utility at a particular time 
t. We also define Viti(qr») = [fa ■ ■ ■ fij] T and Vu 4 (q*) = 
[fii ' ' ' fi.j] T > wnere ai l the vectors are column vectors. 



We define (x) + = max(0,x) and 



The derivative of V along the solution trajectories of the 



x y > 

(x)+ y<0. 



Another way to think of this notation is (x)t = x(l — 
l(-oo,o]( a; )l(-oo,o](2/))> where 1 is the indicator function, i.e.. 
1a(x) = 1 if x € A, and otherwise. 

Motivated by the work in fTT) , we consider the following 
standard primal-dual variable update algorithm: 



iij = k h (fa - Pj)t< > * e J >i e J 



pj 






(8) 



(9) 



Here kf '-, A£ are the constants representing update rates. The 
update rule ensures that, when a variables of interest (q^ or 
Pj) is already zero, it will not become negative even when 
the direction of the update (i.e. quantity in the parenthesis) 
is negative. The tuple (q(t),p(t)) controlled by equations ffify 
and |9} will be referred to as the solution trajectory of the 
differential equations system defined by dSj and (J9l). 

The motivation for the proposed algorithm is quite natural. 
A provider increases its price when the demand is higher than 
its supply and decreases its price when the demand is lower. 
A user decreases his demand when a price is higher than 
his marginal utility and increases it when a price is lower. 
In essence, the algorithm is following the natural direction of 
market forces. 

One key observation is that these updates can be imple- 
mented in a distributed fashion. The users only need to know 
the prices proposed by the providers. The providers only need 
to know the demand of the users for their own resource, and 
not for the resource of other providers (as was the case in 



the analysis of Section IV I. In particular, only user i needs to 
know his own channel offset parameters Cy, j G J . 

The first step to prove the algorithm's convergence is to 
construct a lower-bounded La Selle function V(q(t),p(i)) 
and show that its value is non-increasing for any solution 
trajectory (q(t),p(t)) that satisfies <(8j and |9]). This will 
ensure that (q(t),p(t)) converge to a set of values that keeps 
V(q(t),p(t)) constant. 

B. Convergence of the primal-dual to an invariant set 
We consider the following La Salle function: 



V(q(t),p(t)) = V(t) 



£-~i k q , 



W-qtjW- 



^k] 

3 ■> 



Pj(t) 



W~P*)dP, 



It can be shown that V(q(t),p{t)) > V(q*,p*), i.e., V 
is bounded from below. This ensures that if the function V is 
non-increasing, it will eventually reach a constant value (which 
may or may not be the global minimum V(q*,p*)). 



ov 



system, ^-, denoted by V, is given by: 



dV 






% 



E 



dV 
dp 



-Pj. 



Lemma 7: The value of the La Selle function V is non- 
increasing along the solution trajectory, defined by ([8]) and 
|9}, i.e. V(t) < 0. 



Proof: Proof is given in Appendix C- 1 The proof manip 



ulates the expression for V and shows that it can be reduced 
to the following form: 

i \ 3 , 

+ £(««(*) -««)(/« -!$)• (10) 

i,3 

Using concavity of u^s and properties of the equilibrium 
point q*,p*, we can show that individual elements of the 
summations in (TTOb are non-positive. ■ 

Combining Lemma[7]and the La Salle's invariance principle 
(Theorem 4.4 of fl2) , p. 128) we can prove the following: 

Proposition 8: The pair q(t),p(t) converges to the invari- 
ant set V L = {q{t),p{t) : V(q(t),p(t)) = 0} as t ->• 00. 

It is clear that the invariant set Vl contains the solution 
trajectory that has the value of the unique maximizer of SWO 
(q*(t),p*(t)) = (q*,p*) for all t, since V(q*,p*) = 0. 
However it may contain other points as well. When the 
trajectory (q(t),p(t)) enters the invariant set, it either reaches 
its minimum (i.e., by converging to the unique equilibrium 
point (q*,p*)), or it gets stuck permanently in some limit 
cycle. In either case, the trajectory will be confined to a subset 
of V L - {(q(t),p(t)) : V(q(t),p(t)) = 0}. 

The remainder of this section is to show that the invariant 
set Vl contains only the equilibrium point (q*,p*). This will 
be done in two steps. First, we show that the set Vl has only 
one element for the majority of provider competition instances, 
without any restrictions on the variable update rates. Second, 
we provide a sufficient condition on the update rates so that 
the global convergences to the unique equilibrium point is also 
guaranteed in the remaining instances. 

C. Convergence when providers have decided customers 

In the following two sections we consider the properties of 
the solution trajectory on the invariant set Vl- 

The proof of Lemma [7] shows that individual terms on 
the right-hand side of (JT0|» are non-positive. Combined with 
Proposition [8] we get the following result: 

Corollary 9: On the invariant set Vl, q(t),p(t) are such 
that: 



/ AV3 

3 



(««(*) - ««)(/«(*) - /«) 



=(V Ui ( qi )(t) - V Ui (q*)) T ( qi (t) - q*) = 0,Vz e 2 (11) 
feW - «&)(/« ~Pj) = 0, for all iel,jej, (12) 



where we recall that Vui(qi)(t) — [fa- ■ ■ fij] T and 
Vu i (qt) = [ft 1 ---f* J } T . 

Expressions (fTTb and (TT2j give basic properties of the solu- 
tion trajectories q(t),p(t) on the invariant set. We next prove 
two intermediate results that give further characterization of 
q(t),p{t) on V L . 

Lemma 10: For any point in the invariant set Vl, we have 
fij(t) — f*j for all i E I, j E J . In other words, any user i's 
marginal utility with respect to its demand of any provider j 
equals to the corresponding value at the unique equilibrium. 
In addition, qij(t)(f* i — p*) = on the invariant set Vl- 

The proof is given in Appen dix |C-2 An equivalent way of 
stating the first part of Lemma [10| is that Xi(t) = x* on Vl'. 

^2 1ij(t) c ij = ^2 Qijdj, for a11 * G x - (I 3 ) 



The second part of Lemma 10 claims that q^ (t) (/*, — p* ) = 0. 
From the proof of Lemma ftwe know q*Af*j ~ Pj) = 0- So 
for i,j such that f*. < p*, g* = implies that qij(t) = 0. 
This is good news, since from Lemma [6] most users have zero 
demand to all but one provider at the unique equilibrium. Now 
we know that the same holds on the invariant set Vl. Similarly, 
qij {t) >0 only if /£ =pj. 



Lemma 10 does not preclude the possibility that a demand 
qij(t) for a user with /?• = p* may oscillate between being 
zero and being strictly positive. The following Lemma shows 
that this is not possible: 

Lemma 11: The set Ji(t) = {j E J : qij(t) > 0} does not 
change over time on the invariant set. In addition, pj(t) > 
on the invariant set for all j E J . 

The proof is given in Appendix |C-3| Lemma [TTJ implies that 
if qij(t) = on the invariant set Vl, then <jij(t) = on Vl. 
Also, if q l3 (t) > then q l3 (t) = k^ (/£ ~Pj), as strictly 
positive demand stays strictly positive on the invariant set. 

We are now ready to claim the main result of this subsection: 

Theorem 12: A demand vector of a decided user i E I con- 
verges to the equilibrium demand vector, i.e. lim t _ i . 00 (?,(£) = 
q* . The price pj of any provider j who has at least one decided 
user at the equilibrium (q*,p*), converges to the equilibrium 
price, i.e. lim^oo Pj (t) =p*. 

Proof: Consider an arbitrary user i who at the equilibrium 
has only one preferred provider j, i.e., f~ = p*j for j such that 



q*j > 0, and f*- < p* for j ^ j. By Lemma 10 this implies 



qij(t) > 0, qij (t) = and q^ (t) = for all y^r 3- Combined 



with (13 1, this means that <&j(t) = q% and <jij(t) = 0, i.e., 



user iTdemand vector converges to its equilibrium value. For 
provider j, q^ > so by Lemma [TT] 



=q l3 (t) =k%{f*j- Pj(t))l_ = k%{f*^- p- 3 (t)) 



From this it follows that pj(t) — f* 3 — p|. Further 
differentiation yields pj(t) = 0, meaning that the prices have 
also converged, which completes the proof. ■ 

Theorem 13: If every provider has at least one decided 
customer in the unique equilibrium of the provider competition 
game, the primal dual algorithm converges to this equilibrium. 



Proof: By Theorem 12 the prices of all providers and 



the demand vectors of all decided users converge on the 
invariant set. It remains to be shown that the demand vectors 
of undecided users also converge. By an argument similar to 
the proof of Lemma |2l we can draw a BGR for the undecided 
users. The demand vectors of undecided users that satisfy the 
constraints on the BGR are unique, and thus also converge. ■ 
In most practical cases, where the number of users is 
much larger than the number of providers, all providers will 
have at least one decided user, and hence the convergence 
of the primal-dual algorithm is guaranteed. Next we study 
the more complicated case where some providers do not 
have any associated decided users. In that case, we can still 
prove the convergence of the primal-dual algorithm under mild 
conditions of the variable update rates. 

D. Convergence when providers have no decided customers 

Without loss of generality, we now focus on the problem 
where all the users are undecided at the equilibrium. If we 
can prove that the algorithm converges in this case, then we 
can also prove convergence in the more general case where 
some providers have decided users. Let / be the number of 
undecided users, and J the number of providers. From Lemma 
[6] we know that I < J. 

Theorem 14: Let I < J, and suppose that at the equilibrium 
\{j G J ■ Q*j > 0}| > 1 for all i E 2. The primal-dual 
algorithm converges to the unique equilibrium if the price 
update rates k v , are not integer multiples of each other and the 
demand update rates are equal, i.e. k\- = k Vi G I, Vj E J . 

Proof: We first define some matrices to facilitate the 
proof. Let C be the I x IJ matrix of channel offsets Cy . Let 
I j be the identity matrix of size J x J. Define matrix A to be 
the IJ x J matrix of / vertically stacked identity matrices. Let 
K v = diag(k I ',j E J) be a J x J diagonal matrix containing 
price update rates. Let the IJ x IJ diagonal matrix K q be the 
matrix of demand update rates whose y th entry is KS y = k\ -, 
where i = [(y — 1)/JJ + 1 and j = ((y — 1) mod J) + 1, 
and the off-diagonal entries are all zeros. Here [a^J is the floor 
function representing the largest integer not larger than x. 

By ( p"3"| ), we know that the effective resource of all users 
has converged on the invariant set Vl. We rewrite this as 
cj(t) qi = cfq* for all i € 1, or Cq(t) = Cq*. 

We want to express the primal-dual update algorithm ((HI) 
and Q in matrix form. The final hurdle is getting rid of the 
(x)y operation. From Lemma nil we know that p(t) > and 
that the support sets of vectors q(t) do not change on the 
invariant set. Hence, we can write q^ = fc|- (fy — pj) = 

k ij (fij ~Pj) l(O.oo)fei) (recall that /y = /*, on the invari- 
ant set). This enables us to revise the definition of the update 
rates matrix to be K q — diag(k q A( 0-oo }(qij),i E I,j E J). 
Then, ([8J) and |9]l can be written as: 

q =K" (f - Ap(t)) (14) 

p =KP (A T q(t) - Q) . (15) 

Notice that ([14) and ( fT5J ) form a system of linear equations, 
so the non-linear primal dual dynamics defined in ([8]l and d9l) 
becomes linear on the invariant set. The following result paves 
the way to showing that p(t) is constant on Vl- 



Lemma 15: Let E and B be constant matrices, where B — 

[B; BD- ■■■ ■ BD J - 1 ] (B is a tall matrix), B = CK q A, and 

D = KPA T K<JA. The dimensions of B, B and L> are I J x J, 

JxJ and JxJ, respectively. On the invariant set, Bp(t) = E. 

Proof: The proof is obtained by repeatedly differentiating 



200 



Equations (14i and (15 1 with respect to time. The more 



detailed calculation is given in Appendix C-4 



If we can prove that the rank of B is J, then we could write 
Bjp(t) = E, where Bj is a JxJ matrix constructed by taking 
J linearly independent rows of B. Then p(t) = EBj 1 , which 
implies that p(t) converges on the invariant set. To show that 
the rank of matrix B is J, we use Theorem 6.01 from (JT3J 
(its proof is similar to that of Theorem 6.1 in fT3)): 

Theorem 16: A matrix B = [B; BD; • ■ ■ ; BD^" 1 

column rank if and only if matrix G = „ . _ 

D — Ajii 

column rank for all eigenvalues Xj, j E J of D. 

We now provide a sufficient condition to ensure convergence 

of the primal-dual algorithm. 

Lemma 17: Let matrices B and D be defined as in Lemma 
B 



has full 
has full 



15 



Matrix G 



has full column rank for all 



D-XI 

eigenvalues A of D if k p - ^ ak P , Vj, f E J and all a E N" 
(i.e., as long as the price update rates are not integer multiples 
of each other) and k\, = k, for all i E X,j E J . 



Proof: The proof is given in Appendix C-5 



Combining Lemma 16 and Lemma 17 we see that a unique 
vector pit) can be computed from equation Bp(t) = E, 
meaning that p(t) takes a single value on the invariant set and 
does not change with time. This also means that the demand 
vector q(t) does not change on the invariant set. So since 
p = and q = 0, the primal dual algorithm converges to 
an equilibrium point, which we call a dynamic equilibrium 
point. It can be shown that the dynamic equilibrium point 
is constrained by the same set of equations as the unique 
equilibrium of the provider competition game (refer to the 
proof of Theorem B), Hence, there is only one element in the 
invariant set Vl, and it corresponds to the equilibrium of the 
provider competition game, (q*,p*). This concludes the proof 
of Theorem [14] ■ 

Note that that the condition on the update rates in Theorem 
14 is sufficient but not necessary. In fact, by looking at the 
form of the D matrix from the proof of Lemma 17 we can 



see that a sufficient condition on the price update rates is that 
fc?(53igZ' kij)> where Xj = {i E I : qij > 0}, has a different 
value for each j E J . This condition can be satisfied with 
probability 1, e.g. by drawing fc^'s and fc?-'s independently 
from some continuous distribution. 

VI. Numerical Results and Discussion 

For numerical results, we extend the setup from Example [T] 
where the resource being sold is the fraction of time allocated 
to exclusive use of the providers' frequency band, i.e., Qj = 1 
for j E J . We take the bandwidth of the providers to be 
Wj = 20M Hz, j E J. User i's utility function is a ?; log(l + 
X),=i Qij c ij)i where we compute the spectral efficiency c%j 
from the Shannon formula iW / log(lH — ^r^\hij\ 2 ), qtj is the 
allocated time fraction, E^/Nq is the ratio of transmit power to 
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Fig. 2: Example of equilibrium user-provider association 



thermal noise, and a^ is the individual willingness to pay factor 
taken to be the same across users. The channel gain amplitudes 
\hij\ = -I72 follow Rayleigh fading, where £y is a Rayleigh 

distributed random variable with parameter 1, and a = 3 is 
the outdoor power distance loss. We choose the parameters so 
that the Cy of a user is on average around 3.5Mbps when the 
distance is 50m, and around 60Mbps when the distance is 5m. 
The average signal-to-noise ratio Eb/(NQd a ) at 5m is around 
25dB. We assume perfect modulation and coding choices 
such that the communication rates come from a continuum of 
values. The users are uniformly placed in a 200m by 200m 
area. We want to emphasize that the above parameters are 
chosen for illustrative purposes only. Our theory applies to 
any number of providers, any number of users, any type of 
channel attenuation models, and arbitrary network topologies. 

We first consider a single instantiation with 20 users and 5 
providers. In Fig. [2] we show the user-provider association 
at the equilibrium for a particular realization of channel 
gains, where the thickness of the link indicates the amount 
of resource purchased. The users are labeled by numbers (1- 
20), and the providers are labeled by letters (a-e). This figure 
shows two undecided users (12 and 16), and that certain users 
(1,7,13, and 8) do not purchase any resource at equilibrium. 
For the same realization of channel parameters, Fig. [3] shows 
the channel quality, user demand, and rate at the equilibrium, 
for users 5, 10, 12, 15, and 18. We see that user 15 has 
better channel with provider d than with provider a, but in 
the equilibrium all of his demand is towards provider a. This 
can be explained by looking at the bottom part of Fig. HI 
where the dashed lines indicate equilibrium prices. We see 
that the provider a announces a smaller price than provider d. 
The equilibrium prices reflect the competition among users: in 
Fig. [2] we see that provider b has the most customers, so it is 
not surprising that its price is the highest, as seen in Fig. ffj 

We next consider the convergence time of the discrete time 
version of the primal-dual algorithm. We fix the number of 
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Fig. 3: Channel quality offset, demand and effective 
resource at equilibrium 
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Fig. 4: Evolution of the primal-dual algorithm 
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Fig. 5: Average time to convergence, varying e 



providers to be 5, and change the number of users from 20 
to 100. For each parameter, we run 200 experiments with 
randomly generated user and provider locations and plot the 
average speed of convergence. The convergence is defined as 
the number of iterations after which the difference between 
supply and demand is no larger than eQj. Fig. [5] shows 
the average convergence time for different values of e. In 
general, 200 to 400 iterations are needed for convergence with 
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Fig. 6: Average time to convergence for different 
numbers of providers 
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Fig. 7: Average time to convergence with standard deviation 



e = 10~ 2 , and 100-200 more iterations to get to 10 -3 . In 
Fig. [6] we compare the average convergence time for different 
number of providers. Here we take the stopping criterion to 
be e = 10~ 3 . The convergence time depends on the update 
rates used for the primal-dual: if the update rates are too high, 
then the variables will tend to oscillate, so the algorithm will 
take too long to converge. On the other hand, if the rates are 
too small, the variables may converge too slowly. According 
to our theoretical analysis from Section W\ we randomly 
assigned update variables to ensure global convergence of the 
algorithm. In general a very small or very large number of 
users per provider means that the algorithm will take longer 
to converge. Finally, Fig. I7J shows the average convergence 
time for 5 providers with the standard deviation. The variance 
of convergence time does not change with the number of users, 
except for the case of 20 users. Empirically, the algorithm is 
sensitive to the choice of update rates when the ratio of users 
per provider is smaller than 4. In such cases the prices and 
demands may oscillate and take a long time to converge, which 
can also be seen in Fig. [6] 

A. Discussion 

We first comment on the implication of having undecided 
users at equilibrium. Given a set of prices, the decided users 
can calculate the unique demand vector that maximizes their 
payoffs, while the undecided users have an infinite number 
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of such vectors. In particular, calculating the equilibrium 
maximizing demand vectors for undecided users may require 
cooperation between different providers, which may be chal- 
lenging in practice. 

On the other hand, the number of undecided users is small, 
i.e., not larger than J, and it does not grow with the number 
of users. Future systems may have user action replaced by 
the actions of software agents in charge of connection and 
handover between different providers. In this case, splitting 
over different providers may become feasible. This is similar 
to soft handoff (soft handover), a feature used by CDMA and 
WCDMA standards fl4| . In addition, when the number of 
users is large, the impact of a single user on the price may 
be small. Hence, operating at a non-equilibrium price as the 
result of the decisions of a few undecided users may not have 
a great impact on the experienced quality of service, although 
the exact loss remains to be quantified. 

In the model we consider, we observe the "locally monop- 
olistic" nature of wireless commerce, which does not exist for 
most other traditional goods. Namely, a user that has a strong 
channel to some provider, but a weak one to others, is willing 
to pay a higher price to the provider with the strong channel, 
and is thus not influenced by moderate price changes during 
the price competition. On the other hand, users with similar 
channel gains to all providers will be more sensitive to price 
competition. This local monopoly is in contrast to some other 
wireless resource allocation models where users' association 
is based solely on the price, in which case all users go to the 
provider with the cheapest price (see, e.g. |15|). 

VII. Related Work 

In this paper we have considered a linear-usage pricing 
scheme, which has been widely adopted in the literature (see 
e.g. [4|, |16|). Analyzing such pricing yields various insights: 
for example the existing TCP protocol can be interpreted as 
a usage-based pricing scheme that solves a network utility 
maximization problem |4]. In practice, however, providers 
charge monthly subscription fees. For both voice and data 
plans, these subscriptions are sometimes combined with linear 
pricing beyond a predefined usage threshold. Pure linear pric- 
ing based on instantaneous channel conditions is generally not 
used, although it has received renewed attention due to near- 
saturation of some mobile networks (see, e.g. (17)). Recently, 
AT&T introduced hybrid price plans, consisting of a flat rate 
fee for a certain amount of data, and linear pricing beyond 
that limit (TJ). 

There exists a rich body of related literature on using 
pricing and game theory to study provider resource allocation 
and interactions of service providers. The related research in 
the wireless setting can be divided into several categories: 
optimization-based resource allocation of one provider (e.g., 



T9J-|24|), game theoretic study of interactions between the 



users of one provider (e.g., |25|-|28|), competition of different 
service providers on behalf of the users (e.g., (29), {30]), and 
providers' price competition to attract users (e.g., pT)-(39j). 
Our current work falls into the last category. 

In our work, we have simultaneously considered several 
factors that reflect diverse wireless network scenarios: an 



arbitrary number of wireless providers compete for an arbitrary 
number of atomic users, where the users are heterogenous both 
in channels gains and in willingness to pay. In related work 
where providers price-compete to attract users (31) , [33 1, |34|, 
purchasing a unit of resource from different providers brings 
the same amount of utility to a user; in our work a user's 
utility still depends on the channel gain to the provider. In 
other work where the Wardrop equilibrium concept is used 



(e.g., [15]), users are infinitesimal and non-he terogenous; in 
our work users are atomic and have different willingness to 
pay. One of the early work that explicitly takes into account 
the channel differences for different users on a line is |36|, for 
infinitesimal users and distance-based channel gains. Recently, 
a model similar to ours was used to treat a three-tier system 
(35) , although for specific utility functions. The multiple- 
seller multiple-buyer dynamics in a cognitive radio setting was 
studied in (37) using evolutionary game theory. Finally, |38] 
and 1 39] consider price competition in a multi-hop wireless 
network scenario. 

The design and proof of the decentralized algorithm were 
inspired by Chen et al. (TT) , with several key differences. First, 
their work considers the optimal resource allocation of a single 
OFDM cell. Second, it studies a system where each user has 
a total power constraint. Third, there are often infinitely many 
global optimal solutions in 1 1 1 ] and finally, our convergence 
results are proved with a set of conditions that are less stringent 
than those of (TT) . 

VIII. Conclusions and Future Work 

We provide an overview of the relationship between differ- 
ent concepts used throughout this work in Figure [8] 



Social Welfare Optimization (SWO) 

Section III, Assumptions 1-3 

maximizing vector q* 

Lagrange multipliers p* 

(q*,P*)< 



Provider Competition Game 
Section IV, Assumptions 1-4 

equilibrium price p* 
equilibrium user demand q* 

>(q*,P*) 



Primal-Dual Algorithm, Section V 
lim t _ ) . 0O (qr(t),p(t)) = (<j*,p*) 
>{q*,P*)< 



Fig. 8: Relationship between different concepts 

We have studied the competition of an arbitrary number 
of wireless service providers, who want to serve a group of 
atomic users who are heterogenous in both willingness to pay 
and channel quality. We have modeled this interaction as a 
two-stage wireless provider game, and have characterized its 
unique equilibrium. We have shown that the provider compe- 
tition leads to a unique socially optimal resource allocation 
for a broad class of utility functions and a generic channel 
model. Our results show that some users need to purchase their 
resource from several providers at the equilibrium, although 
the number of such users is upper bounded by the number of 
providers. We have also developed a decentralized algorithm 
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which converges to the equilibrium prices as well as the 
equilibrium demand vectors using only local knowledge. 

Further work may include the study of fractional equilibria 
with the goal of characterizing the losses that occur when 
undecided users are unable to split their resource demand in an 
optimal way. It is also interesting to consider communication 
models where users cause externalities such as interference to 
each other. 

Appendix A 
Social Optimality Proofs 

1. Proof of BGR Property \3\ 

We first examine the properties of the optimal demand 
vector of the SWO problem. We will express the SWO in 
terms of the demand vector q only, by substituting directly 
equation (pi into equation pi. Let p = [p\- ■ ■ pj] be the vector 
of Lagrangian multipliers. The Lagrangian for SWO is then 






3 = 1 V i=l / 



(16) 



It is easy to check that the SWO problem satisfies the Slater's 
condition |40|, and thus the sufficient and necessary KKT 
conditions for an optimal solution (q,p) are as follows: 
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(17) 

(18) 

(19) 

(20) 
(21) 



where with some abuse of notation we use „ Xl ' to denote 

OXi 

dui (x) 

dx _ _sr^J 

X — Xi — 2-,j = l Qij C ij 

The following characterizes the relationship between the 
prices of any two networks from which a user has strictly 
positive demand. 

Recall the support set definition Ji{qi) = {j E J : q%j > 
0}. From ( fP7j ) we see that 



dui(xi) 
dxi 



( 1 8p we can see that 
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> only when - 
mirifcgj ^ is a necessary condition for q^ > for all 
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(22) 



i glj e J . Then, q^ > and gy* > implies - 
minfe G j ^-. In particular, gy > and gy/ > implies 

C ij __ Pj_ 
Ciji Pj' 

We now consider the BGR defined by the support sets 
{Ji}l=i °f undecided users. For any two edges (i, j) and (i, k) 
of BGR, where i is the user index and j, k are the provider 
indices, g,j > and gj/. > so by (22 1 we have ^- = — . 



Suppose that a loop exists in BGR (refer to Fig. [9] for this 
part of the proof). Then, a sequence of nodes ii,ji,i2,J2,---, 
iniJniH exists, where ii,.-.,i n are the user nodes and 
ji,...,j n are the provider nodes, such that (ik>jk) an d 
(ifc-i) ik) are edges in BGR for k — 1, . . .n (with io defined 
as i n ). We assume that the members of the sequence are 
distinct otherwise there is already a smaller loop inside. Since 
both (i k ,jk) and (j k -i,ik) are edges, then ^l 



based on (22 1. A loop in the BGR implies: 



Pk-l 

Pk ' 



Ci^ri Ci 2 \ 
Ciil Ci 2 2 



-m-2 Ci 



PnPl 
PlP2 



Pn-2Pn-l 
Pn-1 Pn 



= 1. 




Fig. 9: A bipartite graph representation loop 



Since 



— is a function of indepen- 



ll Ci 2 2 C i<tl _ in _ 1 

dent continuous random variables, it is also a continuous 
random variable itself. The probability that the product of 
independent continuous random variables equals a constant 
is zero, so we can conclude that a BGR has loops with 
probability zero. In other words, a BGR has no loop with 
probability one. 

2. BGR Algorithm 

Let E be the set of edges, and 1 and J be the set of all 
user and provider nodes, respectively, present in the BGR. The 
demand of undecided users can be found using Algorithm 1. 

Algorithm 1 BGR decoding 

1: For each undecided node i € 1, calculate the checksum 

Pi <- x* 
2: For each provider j £ J calculate the checksum Sj (- 

Qj ~ z2i:(i,j)$g Qij> V? € J 
3: For each q*j > 0, add edge (i,j) to the edge set E 

while E ^ do 

find a leaf node I and associated edge (i,j) 
if the leaf node is a user node then 



% 
else 

Qij 
end if 

Pi <- (Pt 



Si 



and Si 



remove edge (i,j) 
end while 



(Si 



We now give an informal description of an algorithm that 
finds the optimal and unique values of q* for undecided users. 
Since BGR has no loops, it is a (unrooted) tree. Hence, we 
can run a simple iterative algorithm which removes a leaf node 
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(node with a single incoming edge) and its associated edge at 
each iteration. We begin by finding a leaf node in the BGR. 
We then determine the demand of the edge associated to the 
leaf node either from BGR Property (T) or pi. Using this value 
we update the check-sum value of its parent node. Then we 
remove the leaf node and the associated edge. This completes 
one iteration. We repeat the process until there are no more 
edges in the graph. 

The key for Algorithm [T] to work is that the BGR has no 
loops, so a leaf node can always be found in line B] Notice 
that in the last iteration, there will be only one user node i 
and one provider node j left connected by an edge with value 
q* -. The checksums for these two nodes are P, and Sj, which 
satisfy Pi = SjCij since Pi = q*jCij and Sj = q*j. Upon 
completion of the algorithm, the demand of undecided users 
is uniquely defined. 

3. Proof of Theorem H\ 

Assume there exist two optimal demand vectors of SWO 
q* and q'. By Lemma [2] q* and q' have different support 
sets {Si}i=i and {j7/}7=i almost surely. Next, consider a 
convex combination demand vector q x = \q* + Xq' where 
A € (0, 1), A = 1-A. Since x* = £^ =1 qfjCij = £/=i q'ijCij, 

then Y?j=i q ii Ci i = A S/=i Qijdj + A E/=i v'ijdj = <- s ° 
it follows that q x is also a maximizing solution of SWO for 
any A € (0, 1). Then, the support set jMq x ) = {j E J : 

lij = X< £j + ty > °} for user * is $i = 3* U %> for 
all A E (0, 1). In particular, the support sets {J x }\=i are tne 
same for all A E (0, 1), which is a contradiction to Lemma [2] 

Appendix B 
Provider Competition Game Proofs 

1 . Proof of Lemma p] 

It can be verified that UPM satisfies Slater's conditions |40|. 
The necessary and sufficient KKT conditions for an optimal 
solution q { > of UPM of user i are as follows: 



u'iix^dj <pj, j E J 
Qij (u'i(Xi)Cij - Pj) =0, j E J 



j 
where Xi = ^JgijCjj, 



9,>0 



Expression (23 I implies u'^Xi) < a, where a = min^gj- 



(23) 
(24) 

(25) 

Pk 



Based on the utility function of user i, there are two cases: 
u'i(0) < a and 1^(0) > a. 



-Pj < o 



1 



In the first case u£(0)cjj — pj < 0, so u\{Xi 
for all Xi > and for all j E J . This is because the marginal 
utility u^(-) is a strictly decreasing function, by Assumption 
l] Thus, by ( |2"4"| ), q^ — for all j E J . Therefore, qi = and 
by ((25) a;* = 0. So, ((23) -((25) hold for a unique value x* = 0. 

In the second case, u^(0) > a. Then, because u^(-) de- 
creases to zero (Assumption (T), there is a unique Xi > such 
that Ui(^j) = a. We first check that there is a q^ such that 
equations d23)-((25) hold with Xi = X^. Equation d23) holds 



because u'^Xi) = a < Pj/dj for all j E J . Next, by ((24), 
for any j such that Pj/cij > a = u'^xi) we have q^ — 0. For 



any other j, Pj/cij = a = u'^xi), thus, with respect to ( |24) , 
q.jj can take any non-negative value. In particular, for the set 
{j E J : Pj/cij = a} we can choose q^'s so that ( f25) holds. 
Note that gy is positive only when Pj/Cij — a, which 
proves the last part of the lemma. It remains to show that 
Xi is the only value of Xi for which a q { satisfying ((23)- 
((25) exists. For any Xi < Xi, u'^Xi) > a so ((23) is violated 
for j E aigmin pk/cik- For any Xi > Xi, Ui{xi) < a so 
Ui(Xi)Cij — Pj < for all j E J . Then, ( |24) implies that 
qij = for all j E J , meaning that Xi — 0, which contradicts 
Xi > Xi > 0. Therefore Xi is the unique searched value x*. 

2. Proof of Theorem [5] 

Assume that the providers charge prices p = [pi . . . p j] to 
the users. Then, each user faces a local maximization problem 
UPMj(p), as defined in ((6). 

Equations d23) - d25), together with q { > for all i El, are 
equivalent to equations (fl7)-(fT9) and (21 1, Furthermore, under 
assumption HJ for q to be an SPE of the provider competition 
game, the demand to each provider must equal its supply, i.e., 
2~2i=i Qij — Qj f° r a ^ 3 S J (i-e., equations (|20)). Hence, the 
SPE is a price vector tuple (p, q) that satisfies KKT conditions 
(T7)-((2T). But, the KKT conditions ((T7)-((2T) are necessary and 
sufficient for any vector tuple (p, q) to be the maximizing 
solution of SWO. Hence, we have established formal equiva- 
lence between the SPE of the provider competition game and 
the maximizing demand vector and Lagrangian multipliers of 
the SWO problem (q* ,p*). Hence, (p* , q*) forms the unique 
SPE of the provider competition game. 

3. Proof of Lemma [6] 

Given an arbitrary price vector p, we can construct the 
corresponding BGR(p) from the users' preference sets. First 
consider a BGR that is connected (single component). We start 
drawing the graph with a single undecided user node, and then 
add the provider nodes that are connected to this user node. 
There should be at least two such provider nodes. Then we 
add another undecided user node which shares one common 
provider node with the existing undecided user node. This new 
undecided user node will bring at least one new provider node 
into the graph, otherwise it leads to a loop in the graph. We 
repeat this process iteratively. Since the number of provider 
nodes is upper-bounded by J, the total number of undecided 
user nodes is upper-bounded by J — 1. The "— 1" is due to 
the fact that the first undecided user node is connected to (at 
least) two new provider nodes in the graph. If we consider a 
BGR with b disconnected subgraphs, it can be shown that the 
total number of undecided users is bounded by J — b. 

Appendix C 
Primal Dual algorithm supporting proofs 

1. Proof of LemmaVA — Proof that V < 0: 

For the optimal demand vector q* of the SWO problem and 
the associated Lagrange multipliers p* , we see that /*■ = p* 
whenever q*. > from equation (fTS), and /* < p* when 
q*j = from equations ( 17 1 and (II 8b. Similarly, equation ( 18 1 



13 



ensures that f*- < p* implies g* = 0. This fact will be used 
shortly. For our La Salle function: 



V = 



(a) 



- £(«« - Q*j)(fij - Pi)tn + X^ ~ Pj)(%2 9a - Qj) 



p.i 



(b) 



< X]fej ~ 9*j)(fij -pj) + ^(pj - Pj)(52 *■? ~ QA 



»>j 



(c) 



X)(X^ ~ ««)(/« - ft])) + Xfe ~ 9*3) (/t* 



!-J 



where (a) follows from the definition of V and V, (b) can be 
readily verified by examining all the cases, and (c) is obtained 
by some algebraic manipulation. The expression for V is now 
in such a form that we can prove V < 0. 

First, the following is true for any two vectors q x and q 2 
due to concavity (see, e.g. Section 3.1.3 of pO)): 

Vuf (q x )(q 2 ~ Qi) > Ui{q 2 ) - ^(<?i), (26) 

Vwf (g 2 )(«i - Q2) > Ui(li) ~ Ui{Qi)- ( 27 ) 

Substituting q x = q { (t) and q 2 = q* in (|26j and (|27j gives: 
Vuf( qi )(q* - qA > Vu*( qi )(q* - qA, 
which can be rewritten as 



3 



(?«(*) " <&)(/«(*) - /«) < 0. for a11 * € *> (28) 



which we recognize as one of the components of the first 
term in ( fl~0| >. Considering a component of the second term in 

([To), (gy(t) - gfc)(/« - Pp> we first reca11 from the KKT 

conditions ( fl~8] > that either /*■ = p*, in which case (tfoj — 
qtj)(fij ~Pj) = 0, or f* < p*, in which case q*j = so 
(?« ~ Qij)(f*j ~P*j)< 0. Hence < 

( ffy (<) - <?*• )(/*•- P*) < 0, for all i e X, j e J". (29) 

which completes the proof that V < 0. 



The second part of the Lemma is simpler to prove. Equation 

(12]) gives g« (*)(/£-$) = «&(/« _ Pp- which is ec l ual t0 
zero by equation ([18). This concludes the proof. 



5. Proof of Lemma 11 — Proof that the support set of primal 
and dual variables does not change on the invariant set 

We prove that the set of positive demands does not change 
over time on the invariant set by contradiction. Suppose that 
qtj (t) > but qij (t + r) = for all r such that < r < e, 
where e is a small number. Then qijit) < but q i3 (t + r) = 0, 
so we have 

lim V] qij(t)cij - qij(t + r)cij < 

On the other hand differentiating equation ( p"3[ ) with respect to 

time yields Y.^j, 9«(*) c y = ° and Eje^ 1ij(t+T~)cij = 0, 
which is a contradiction. So a non-zero g.y stays non-zero. 
Suppose now that qijit) = but %j(i + Tif) > for some 
Tij > 0. After time ry variable g^- becomes non-zero and 
stays that way forever (according to the argument we made 
earlier in the proof). Then no q^ escapes from the boundary 
after time t + T*, where r* = max^j ry. 

Now we prove the second part of the Lemma. Similar to 
the previous argument, once the pj (t) > on the invariant set, 
Pj (t + r) > for all r > 0. It remains to show that pj (t) > 
on the invariant set. But, if pj (t) = were true, then q^ > 
would imply lim t _ i . 00 qij(t) — 00, which would violate ( 13 1. 



4. Proof of Lemma [75] — Proof that Bp = E 

Since Cq(t) — Cq* is a constant on Vl, then taking the 
time derivative of both sides gives Cqit) = 0. Then: 

Q = Cq(t) = CK c '(f*-Ap(t)) 

^ck q Ap{t) = ck q f\ 

which means that CK q Apit) equals a constant on the invari- 
ant set. Taking the time derivative one more time yields 



2. Proof of Lemma 10 — Proof that marginal utility on Vl is 
the equilibrium marginal utility. 

Recall that for user i, fij{qA — U q^ Ci 



all j e J, where Xi{t) — YljQij(t) c ij an d u 'i( x i(t)) i s a 
scalar function of time. Hence, 



= CK q Ap(t) =CK q AK p (A T q(t) - Q) 
■& CK q AK p A T q{t) =CK q AK p Q = constant 

Repeating the derivative operation 2n + 1 times will yield 
u5(a?i)cy for ck q A(K p A T k q Ay i p=Ck q A{K p A T k q A) n - 1 K p A T k q f, 
= CK q AD n p =CK q AD n - 1 K p A T k q f* = const 



V«f(qfi) = u-(xi)[cii • • • Cij] =u' i {x i )cj. 

We now consider equations (fTTh, which can be rewritten as: 

Vuf (g;)(<Zj - q^ =Vuf(q*)( qi - q*), or, 
Ui(xi)ci(qi - q*) =u'i{x*)cf(qi - q* ), 

which leads to u'^fx^Axiif) — u' i (x*)Axi(t), where 
Axi(t) = cf(q t (t) - q*). Then, either Ax t (t) = 0, or 
Ui(Xi(t)) = u'^x*). The necessary and sufficient condition for 
both is that that Xi(t) = x*, for q.^t), ielon the invariant 
set Vl- An immediate corollary is that fij(t) = f*j on Vl for 
all iel,j ej. 



where we defined D = K p A T K q A. Note that D n here stands 
for "D to the n th power". Let B = CK q A, then we can write: 

[B BD ■■■ BD 1 - 1 ] T pit) =Bp(t) = constant, 
which completes the proof. 



5. Proof of Lemma 1 7 



Consider matrix D = K p A T K q A. It can be verified that 



D 







^(E, 



ex, 



k 



u 
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(a) 



*?|Zl| 











kl\l 2 \ 





k P j\Tj\ 



where X 3 ■ = {i G X : q^ (t) > 0} is the set of users of provider 

k. Without loss of generality, 



j, where (a) follows from kh 



take k — 1. Now, assuming that fc^'s are not integer multiples 
of each other, we can see that at most one row of D — XI can 
be all-zero vector for any eigenvalue A of D (indeed, exactly 
one since D is a diagonal matrix). So, we can always choose 
J — 1 rows of D — XI that are non-zero and also linearly 
independent. An example of G' is: 



fc?|Zi|-Aj 







*S 



a) 



K i.J- 




^{Xj-il-Xj 

CiJ-l^M+{0.iJ- 











) k h c iJ 



Without loss of generality, assume that we choose rows 1 
through J — 1 (i.e., we consider the J th eigenvalue Xj = 
k P j\Ij\). Now, consider matrix B 1 — CK q A. It can be verified 



Hi 



1,3 

,1 



Sj^hJ- 



ij^{0,oo)\Hij 



(qij(t)). Since qij(t) > for at least one i = i, 



'jCij l(o, oo) (qi.J (*)) = k?jCij, so the i row has an entry at 
the J th column. We then append this row to the J— 1 rows we 
took from D — Ajl, forming a submatrix of G, denoted G". 
The submatrix G" is then a lower triangular matrix. Lower 
triangular matrices have full rank, in this case J. Since the 
submatrix of G has the full column rank, then G also has the 
full column rank J. The procedure works in the same way for 
any eigenvalue Xj of D we may choose, since we can always 
relabel the providers. This completes the proof. 



References 



[i] 



[21 



[3] 



"FCC press release: FCC adopts rules for unlicenced use of television 

white spaces," November 2008. 

S. Leyffer and T. Munson, "Solving multi-leader-follower games," 

Preprint ANUMCS-P 1243 -0405, 04 2005/04/22. 

J. Wang, D. Chiu, and J. Lui, "A game-theoretic analysis of the 

implications of overlay network traffic on ISP peering," Computer 

Networks, vol. 52, no. 15, pp. 2961-2974, October 2008. 
[4] F. Kelly, A. Maulloo, and D. Tan, "Rate control for communication 

networks: Shadow prices, proportional fairness and stability," Journal 

of the Operational Research Society, vol. 49, pp. 237-252, 1998. 
[5] X. Lin and N. B. Shroff, "Utility maximization for communication 

networks with multipath routing," IEEE Transactions on Automatic 

Control, vol. 51, no. 5, pp. 766-781, 2006. 
[6] T. Voice, "Stability of congestion control algorithms with multi-path 

routing and linear stochastic modelling of congestion control," Ph.D. 

dissertation, University of Cambridge, Cambridge, UK, May 2006. 
[7] D. P. Bertsekas, Nonlinear Programming. Athena Scientific, Sept. 1999. 
[8] D. Fudenberg and J. Tirole, Game Theory. MIT Press, 1991. 
[9] A. Mas-Colell, M. D. Whinston, and J. R. Green, Microconomic Theory. 

Oxford University Press, 1995. 
[10] A. Rad, J. Huang, M. Chiang, and V. Wong, "Utility-optimal random 

access without message passing," IEEE Transactions on Wireless Com- 
munications, vol. 8, no. 3, pp. 1073-1079, 2009. 
[11] M. Chen and J. Huang, "Optimal resource allocation for OFDM uplink 

communication: A primal-dual approach," in CISS, 2008, pp. 926-931. 
[12] H. Khalil, Nonlinear Systems. Prentice Hall, 2002. 
[13] C.-T. Chen, Linear System Theory and Design, 3rd ed. Oxford 

University Press, 1998. 
[14] D. Wong and T. J. Lim;, "Soft handoffs in CDMA mobile systems," 

IEEE Personal Communications, vol. 4, no. 6, pp. 6-17, December 1997. 



[15 
[16 
[17 

[18 

[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 

[27 

[28 
[29 

[30 
[31 

[32 

[33 

[34 

[35 
[36 
[37 

[38 
[39 
[40 



P. Maille and B. Tuffin, "Analysis of price competition in a slotted 
resource allocation game," in INFOCOM, 2008, pp. 888-896. 
H. Shen and T. Basar, "Differentiated internet pricing using a hierarchi- 
cal network game model," ACC, vol. 3, pp. 2322-2327, lanuary 2004. 
Fierce Wireless, "Is usage based pricing in- 

evitable?" [Online]. Available: http://www.fiercewireless.com/story/ 
usage- based- pricing- inevitable/20 1 0-02- 03 

Engadget, "AT&T makes sweeping changes to 

data plans, iPhone tethering coming at OS 4 
launch." [Online]. Available: http://www.engadget.com/2010/06/02/ 
atandt-makes-sweeping-changes-to-data-plans-iphone-tethering-comi/ 

C. U. Saraydar, N. B. Mandayam, and D. J. Goodman, "Efficient power 
control via princing in wireless data networks," IEEE Transactions on 
Communications, vol. 50, no. 2, pp. 291-303, February 2002. 

P. Marbach and R. Berry, "Downlink resource allocation and pricing for 
wireless networks," in INFOCOM, vol. vol. 3, 2002, pp. 1470-1479. 
M. Chiang and J. Bell, "Balancing supply and demand of bandwidth in 
wireless cellular networks: Utility maximization over powers and rates," 
in INFOCOM, vol. 4, 2004, pp. 2800-2811. 

D. Acemoglu, A. Ozdaglar, and R. Srikant, "The marginal user principle 
for resource allocation in wireless networks," in CDC, vol. 2, 2004, pp. 
1544-1549. 

J. Musacchio and J. Walrand, "WiFi access point pricing as a dynamic 
game," IEEE/ACM Transactions on Networking, vol. 14, no. 2, 04 
2006/04/01. 

L. Jiang, S. Parekh, and J. Walrand, "Base station association game in 
multi-cell wireless networks," WCNC, pp. 1616-1621, 2008. 
S. Adlakha, R. Johari, and A. J. Goldsmith, "Competition in wireless 
systems via Bayesian interference games," CoRR, 2007, informal pub- 
lication. 

R. Etkin, A. Parekh, and D. Tse, "Spectrum sharing for unlicensed 
bands," in DySPAN, 2005, pp. 251-258. 

J. Huang, R. A. Berry, and M. L. Honig, "Distributed interference 
compensation for wireless networks," IEEE Journal on Selected Areas 
in Communications, vol. 24, no. 5, pp. 1074-1084, 2006. 
D. Kalathil and R. Jain, "Spectrum sharing through contracts," in 
DySPAN, 2010. 

C. Zhou, M. L. Honig, and S. Jordan, "Utility-based power control 
for a two-cell CDMA data network," IEEE Transactions on Wireless 
Communications, vol. 4, no. 6, pp. 2764-2776, 2005. 

L. Grokop and D. N. Tse, "Spectrum sharing between wireless net- 
works," in INFOCOM, April 2008, pp. 201 - 205. 
A. Zemlianov and G de Veciana, "Cooperation and decision-making 
in a wireless multi-provider setting," in INFOCOM, vol. 1, 2005, pp. 
386-397. 

S. Sengupta, M. Chatterjee, and S. Ganguly, "An economic framework 
for spectrum allocation and service pricing with competitive wireless 
service providers," in DySPAN, 2007, pp. 89-98. 

J. Jia and Q. Zhang, "Competitions and dynamics of duopoly wireless 
service providers in dynamic spectrum market," in MohiHoc, New York, 
NY, USA, 2008, pp. 313-322. 

O. Ileri, D. Samardzija, T. Sizer, and N. B. Mandayam, "Demand 
responsive pricing and competitive spectrum allocation via a spectrum 
server," in DySPAN, 2005, pp. 194-202. 

J. Acharya and R. D. Yates, "Service provider competition and pricing 
for dynamic spectrum allocation," in GameNets, 2009, pp. 190-198. 
H. Inaltekin, T. Wexler, and S. B. Wicker, "A duopoly pricing game for 
wireless IP services," in SECON, 2007, pp. 600-609. 

D. Niyato, E. Hossain, and Z. Han, "Dynamics of multiple-seller and 
multiple-buyer spectrum trading in cognitive radio networks: A game- 
theoretic modeling approach," IEEE Transactions on Mobile Computing, 
vol. 8, no. 8, pp. 1009-1022, 2009. 

Y Xi and E. M. Yeh, "Pricing, competition, and routing for selfish and 

strategic nodes in multi-hop relay networks," INFOCOM, pp. 1463- 

1471, 13-18 April 2008. 

L. Gao and X. Wang, "A game approach for multi-channel allocation in 

multi-hop wireless networks," in MobiHoc '08. New York, NY, USA: 

ACM, 2008, pp. 303-312. 

S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge 

University Press, 2004. 



